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Isothermal profiles of the extended meniscus in a quartz cuvette were measured in the 
earth ’s gravitational field using an image-analyzing interferometer that is based on com- 
puter-enhanced video microscopy of the naturally occurring interference fringes. These 
profiles are a function of the stress field. Experimentally, the augmented Young- Laplace 
equation is an excellent model for the force field at the solid -liquid -vapor interfaces 
for heptane and pentane menisci on quartz and tetradecane on SFL6. Effects of refrac- 
tive indices of the solid and liquid on the measurement techniques were demonstrated. 
Experimentally obtained values of the disjoining pressure and dispersion constants were 
compared to those predicted from the Dzyaloshinskii- Lifshitz- Pitaevskii theory for an 
ideal surface and reasonable agreements were obtained. A parameter introduced gives a 
quantitative measurement of the closeness of the system to equilibrium. The nonequtlib- 
rium behavior of this parameter is also presented. 


Introduction 

In thin liquid films, the shape-dependent interfacial stress 
field controls fluid flow and heat transfer (Derjaguin et al., 
1965; Huh and Scriven, 1971; Potash and Wayner, 1972; 
Miller and Ruckenstein, 1974; Wayner et al„ 1976; Teletzke 
et al., 1987; de Gennes, 1985; and Wayner, 1991). These con- 
cepts are particularly useful under microgravity conditions 
where the interfacial force field can predominate in systems 
much larger than those in the earth’s gravitational field. When 
the gravitational body force is essentially removed, the shape 
of a constrained liquid volume with a free interface in a con- 
tainer changes dramatically to reflect the new force field. The 
resulting equilibrium shape under microgravity conditions de- 
pends on the intermolecular force field that changes rapidly 
in the vicinity of the liquid-vapor and liquid-solid interfaces. 
This liquid film shape is controlled by fixing the container 
shape, the liquid and solid substrate surface properties, and 
the volume of the liquid. Under nonequilibrium conditions, 
microgravity fluid dynamics and change of phase heat trans- 
fer are a function of further changes in the shape of the fluid 
volume from its equilibrium shape. Herein, we discuss the 
results of ground-based equilibrium experiments that are a 
precursor to nonequilibrium microgravity experiments. We 
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note that the experimental cell is also of general use in deter- 
mining intcrfacial forces. 

In particular, we are concerned with the experimental study 
of the generic constrained vapor bubble thermosyphon 
(CVBT) system presented in Figure 1 for a microgravity envi- 
ronment. For a completely wetting system, the liquid will coat 
all the walls of the chamber. For a finite contact angle sys- 
tem, some of the walls will have only a small amount of ad- 
sorbed vapor that changes the surface properties at the 
solid-vapor interface. Liquid will fill a portion of the comers 
in both cases. If temperature at End (2) is higher than End 
(1), because of an external heat source and sink, energy flows 
from End (2) to End (1) by conduction in the walls and by an 
evaporation, vapor flow, and condensation mechanism. The 
condensate flows from End (1) to End (2) because of the in- 
termolecular force field that is a function of the film profile. 
A CVBT with a very small cross section has been called a 
micro heat pipe (Cotter, 1984; Peterson, 1992). The promise 
for high heat flux cooling has generated renewed interest in 
the application of micro heat pipes as efficient cooling de- 
vices. Potential applications are many, namely, cooling of 
electronic circuitry and thermal control of spacecraft. From a 
more fundamental point of view, the CVBT conveniently per- 
mits the study of interfacial transport concepts. 

There is a “pressure jump’’ at the liquid-vapor interface, 
due to the anisotropic stress tensor near interfaces. For many 
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perimentai results of a study of the three-dimensional ex- 
tended meniscus in a configuration of general importance to 
the microgravity environment: the constrained vapor bubble 
thermosyphon. Experimentally, we find that the augmented 
Young-Laplace equation is an excellent model for the force 
field at the solid-liquid-vapor interfaces for heptane and 
pentane menisci on quartz and tetradecane on SFL6. 

The effects of both the liquid-vapor and liquid-solid in- 
terfaces on the effective pressure jump at the liquid-vapor 
interface of the extended two-dimensional meniscus have 
been modeled using the following augmented Young-Laplace 
equation: 


B B 
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years, the classic Young-Laplace equation of capillarity has 
been successfully used to describe the pressure jump at a 
curved liquid-vapor interface. An example of its use to de- 
scribe the fluid dynamics in a micro heat pipe is given in a 
publication by Wu and Peterson (1991). In this case, the pres- 
sure jump is a function of the liquid-vapor surface tension 
and the interfacial radius of curvature. However, near the 
liquid-solid interface, additional changes in the stress field 
within the liquid occur because of changes in the intermolec- 
I.iar force field due to solid molecules replacing liquid 
molecules. This leads to the augmented Young-Laplace 
model for the pressure jump at the liquid-vapor interface 
(Deijaguin and Kussakov, 1939; Potash and Wayner, 1972; 
Blake, 1975; Derjaguin and Churaev, 1976; Teletzke et al., 
1976; Renk et al., 1978; Kamotani, 1978; Holm and Goplen, 
1979; Moosman and Homsy, 1980; dc Gennes, 1985; Truong 
and Wayner, 1987; DasGupta et al., 1993a). These long-range 
van der Waals forces have been found to be extremely impor- 
tant in that they lead to the concept of an extended evaporat- 
ing meniscus (Potash and Wayner, 1972; Wayner et al., 1976). 
In a completely wetting evaporating system, a thin adsorbed 
film extends for a long distance beyond the classic equilib- 
rium meniscus (such as Derjaguin et al., 1965; Wayner, 1991). 
The thin film controls the important processes of spreading 
and wetting. It forms a thin liquid bridge between the classi- 
cal menisci ’ formed in the comers of the chamber presented 
in Figure 1. Recently, Khrustalev and Faghri (1993) and 
Swanson and Peterson (1993) have used this type of model to 
analyze a micro heat pipe. Heat transfer from a stable evapo- 
rating thin film in the neighborhood of a contact line was 
analyzed by Brown et al. (1993). Swanson and Herdt (1992) 
have used the three-dimensional augmented Young-Laplace 
equation to develop a mathematical model describing the 
evaporating meniscus in a capillary tube. On the other hand, 
experimental results have been presented that question the 
validity of the augmented Young-Laplace equation: Using 
aqueous films, Lam and Schechter (1991) found deviations 
from the predicted profile obtained using a technique pre- 
sented by Hirasaki (1990). Herein, we present the initial ex- 


where P, is the liquid pressure, P u is the vapor pressure, A is 
the Hamaker constant (negative for a completely wetting liq- 
uid), B is the retarded dispersion constant for thicker films, 
8(x) is the film thickness, K is the curvature, and <r is the 
surface tension. The first term on the righthand side of Eq. 1 
is called the disjoining pressure, - II, and it represents the 
change in the body force on the liquid due to the long-range 
van der Waals forces between the liquid and solid over a nar- 
row range of thicknesses. 

According to Eq. 1, the effective pressure in the liquid is 
reduced below that in the vapor by both capillarity and dis- 
joining pressure in a completely wetting system. This leads to 
a reduction in the vapor pressure. However, the vapor pres- 
sure reduction can be offset by a temperature increase to 
obtain the vapor bubble thermosyphon presented in Figure 1. 
The chemical potential field is a function of the temperature 
and the effective liquid pressure that is a function of the shape 
of the liquid film. Therefore, the process is understood by 
measuring the temperature field and the liquid film shape. 


Experiments 

A diagram of the cell and the experimental setup is shown 
in Figure 2. The readily available UV fused quartz cell (Type 
607 Fluorimeter cell, NSG Precision Cells, Inc.) was essen- 
tially a small cuvette, selected to facilitate optical observation 
and measurement of the liquid meniscus. The cell was square 
(3 mm x 3 mm inside) in cross section and about 40 mm 
long. The SFL6 glass cell was specially made from optically 
flat glass and is not readily available. The cell was cleaned 
after manufacturing and delivered in a closed container. The 
cells were claimed to be ready for spectroscopy. The cell was 
opened only inside a class- 100 clean hood, rinsed repeatedly 
with pure ( + 99%, as supplied by Aldrich Chemicals) and fil- 
tered test liquid before partially filling it. We examined the 
inside surface of the cell using high power microscopes and 
found them uniformly clean and sufficiently free of dust par- 
ticles. As a result, we decided to use the quartz cell without 
any further wet cleaning. At the end of each experiment the 
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Figure 2. Experimental setup. 

The source of illumination was through the objective of the 
microscope. 


cell was heated inside the clean hood at about 300°C for an 
hour and then rinsed with the new test liquid. We note that 
slight dehydration of the surface may take place during heat- 
ing, but since we used only nonpolar liquids (alkanes), that 
would not affect adsorption on quartz (Gee et al., 1989). Us- 
ing this procedure we were able to obtain a satisfactorily clean 
surface for experiments. The present study was conducted as 
a precursor to the nonisothermal CVBT experiments, mainly 
to develop the method for accurately measuring the liquid 
film thickness profile and to estimate the dispersion constant, 
in situ , at the start of the experiment. The isothermal results 
characterize the interfacial force field and the information is 
critical to the operation and the heat sink capability of the 
CVBT. The experimental cell was kept inclined to demon- 
strate the changes in shape of the meniscus as a result of 
externally imposed conditions (gravitational field, in this 
case). Examples of two cross sections of the cell that are a 
function of the system, location, gravitational field, volume of 
liquid, and cell size are presented in Figure 3. 

The film thickness profile near the comers of the experi- 
mental cell was measured using image analyzing interferome- 
try (IAI). A drawing of this section of the cell is presented in 
Figure 4. Interference phenomena were used to determine 
the profile of the capillary meniscus in the thickness range 
8 £ 0.1 yn n. A detailed description of the IAI techniques and 
hardware for film thickness profile measurement were pre- 
sented in Sujanani (1990), Sujanani and Wayner (1991), and 
DasGupta et al. (1991,1993a,b) for a reflecting substrate. 
Blake (1975) used a transparent substrate in a related study. 


Figure 3. Example of tho effect of gravity on the menis- 
cus shape. 


Figure 4. Location of microscope. 

The image captured from the microscope through a CCD 
camera was digitized into 490 (vertical) x512 (horizontal) 
space pixels and assigned one of 256 possible gray values rep- 
resenting an intensity from 0 (black) to 255 (white). Thus each 
microscopic pixel acted as individual, simultaneous light sen- 
sors. For the present experimental setup, each pixel repre- 
sented a diameter of 0.3 /am. 

The cell was partially filled up with the liquid and placed 
on the microscope stage. The whole assembly was tilted with 
respect to the horizontal at two different angles in the two 
settings used in this study. The system was left to equilibrate 
with the surroundings for 3 h before taking any data. An op- 
tical interference pattern representing the thickness profile 
was readily observed (Figure 5). The pictures of the interfer- 
ence patterns formed at the comers of the cell at two specific 
points (A and B, as shown in Figure 2) for each angle setting 
were captured in the image processor, The exact locations of 
A and B varied with each system so that interference from 
small particles adhering to the surface, which caused nonpar- 
allel fringes in some regions, could be avoided. The vertical 
separation in the gravitational field between locations A and 
B was Ay, which is given below. From each image a plot of 
the pixel gray value vs. distance was extracted. The gray value 
at each pixel was a measure of the reflectivity. As is evident 
from Figure 5, the reflectivity underwent a cyclic change with 
increase in film thickness. The computer program scanned 
the peaks and valleys and filtered the noise from the real 
peaks/valleys. It then interpolated peak/valley envelopes and 
by analyzing the relative reflectivity of any pixel with respect 
to these (dark and light pixel envelopes) determined a film 
thickness at every pixel (DasGupta et al., 1991), The fact that 
the extended capillary meniscus merged smoothly to an ad- 
sorbed flat film was utilized to estimate the adsorbed film 
thickness from the gray value data and the peak/valley en- 
velopes. The adsorbed film is obvious in the gray value plot 
since G(jc)<G max (x) in the region x > 50 pm. An estima- 
tion of the overall error associated with the measurements 
and data reduction procedures can be obtained by viewing 
the data in the comparison between theory and experiments 
presented below in Figure 11. We note that independent 
measurements of the curvature in the thick film region and 
the adsorbed film thickness in the flat region are needed. 

Figures 6 and 7 are two examples of the measured film 
thickness profiles with heptane as the test liquid. Each pic- 
ture was for a specific angle of inclination. In an isothermal 
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Figure 5. Interference fringes and the gray value plot. 


system of spreading liquid in an inclined cuvette, the curva- 
ture at a constant height in the gravitational field should re- 
main constant in the region where dispersion forces can be 
neglected. So the film profile in this range approximated a 
parabola and a plot of 5 1/2 vs. x was a straight line. As can 
be seen from Figures 6 and 7, 8 1/2 vs. x was nearly a straight 
line in all the cases, showing the proximity of the cases to 
isothermality. On the other hand, it is clear from these fig- 
ures that the thickness profiles at point B in each of these 
angle settings was steeper than at point A as a result of higher 
hydrostatic pressures and subsequent higher curvatures. 
Again a comparison of Figures 6 and 7 demonstrates that, as 
expected, an increase in angle resulted in an increase in the 
film curvature. Plotting 5 1/2 vs. x , therefore, clearly illus- 
trates the sensitivity of the film thickness shape to externally 
imposed conditions. This will not be obvious from a plot of 5 
vs. x. In addition, the adsorbed film thickness, S Q , is clearly 
seen in the region x £ 50 fim. 

Refractive index and contrast 

The contrast of the interference fringes of a liquid film 
formed on a solid substrate is a strong function of the refrac- 
tive indices of the two media, more specifically, of the differ- 
ence of the refractive indices of the solid and the liquid. Since 
we used optical techniques (microcomputer-enhanced video 
microscopy of the interference fringes) for the measurement 
of the film thickness profiles, the relative contrast of dark 
and bright fringes has a profound effect on the measurement 
capability and its accuracy. If the two refractive indices are 


FILM THICKNESS PROFILES 

Liquid : Heptane 
Angle of Inclination : 6.24 



50 100 

DISTANCE ( |im) 


Figure 6. Film thickness profiles of Heptane at different 
positions (point A and B) for an angle of incli- 
nation equal to 6.24°. 


close to each other, it will be difficult to see the fringes and 
perform any image analysis due to the high noise in the cap- 
tured picture. The following discussion gives a qualitative 
treatment of the problem. 

The relative contrast, C, is defined in the following way, 


* l +*2 2 

r,-r 2 • 


( 2 ) 


Here R } and R 2 are the Fresnel coefficients that depend on 
the refractive indices of the solid, liquid and the vapor: 



Figure 7. Film thickness profiles of Heptane at different 
positions (point A and B) for an angle of Incli- 
nation equal to 9.92°. 
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Figure 8. Relative contrast vs. liquid and solid refrac- 
tive indices. 


K, 


n s -n. 


*2 




n { + n 0 


(3) 


Figure 8 is a plot of the relative contrast as a function of the 
solid and liquid refractive indices. As is evident from the fig- 
ure, a drastic reduction in contrast takes place when the re- 
fractive index of a liquid being tested in the quartz cuvette is 
dose to the refractive index of quartz (n s - 1.5). We observed 
noticeable reduction in sharpness and contrast of the fringes 
between pentane («/ - 1.358) and heptane (n, - 1.387) 
menisci on quartz. But for tetradecane (n { - 1.429) in our 
quartz CVBT, the fringes were blurred and the noise level 
was unacceptably high even after using image-enhancement 
techniques. Numerical values of n, R y and C are given in 
Table 1. 

To solve this problem, we designed another CVBT with 
the same dimension but using a spedal high refractive index 
glass (SFL6, n s - 1.8126 at 548-nm wavelength) to enhance 
the contrast and ensure the accuracy and ease of the optical 
measurements. We observed dramatic improvements in con- 
trast (as can be expected from the results presented in Figure 
8) and were able to conduct experiments with higher refrac- 
tive index liquids (such as tetradecane). 


where both effects are present. For the isothermal none\apo- 
rating (additional heat input, Q = 0) cases considered here 
the liquid pressure will remain the same, irrespective of the 
position at a constant gravitational level. In our experiments, 
the following equation applies at a fixed hydrostatic head: 

B 

crK 7 = <jK^ (? = 0. (4) 

<5 4 

We note that the curvature at the thicker portion of the 
meniscus (A x ) is nearly constant in the region studied 
(DasGupta et aL, 1993b). In Eq. 4, the curvature (K) can be 
expressed in the following form. 


d 2 S 



(5) 


Using a simplified form of curvature, which is valid only if 
the square of the slope is small compared to one, Eq. 4 can 
be written in the following way: 


d 2 S B 

,r d^~¥~ <TK ^ 


( 6 ) 


For the isothermal experiments reported in this study the 
value of the maximum slope is 0.15, and hence using the sim- 
plified form of curvature is justified. The following variables 
are introduced next to modify Eq. 6 to obtain Eq. 8: 


n- 


8 



dS -B \ 1 
dZ 1 + ( cr K x 5o ) v* 


1. 


(7) 

( 8 ) 


Equations 4-8 are valid for an equilibrium situation where no 
evaporation or condensation is taking place. A dimensionless 
variable* a, is defined next. 


- B 

4 — 

ctK^Sq ' 


(9) 


Theory 

The augmented Young- Laplace equation can be written 
for a point in the thicker portion of the meniscus, where the 
disjoining pressure effects are negligible and for another point 


Equation 8 can now be written as 


dZ 2 


( 10 ) 


Table I. Relative Contrast for Various Liquid-Solid Systems 


Liquid 



R 2 

Contrast (C) 

Quartz SFL6 
n,~l . 5 *,-1.813 

Quartz 
/t T * 1.5 

SFL6 
n s -1.813 

Pentane 

1.358 

0.0497 0.1435 

0.1518 

3.895 

1.255E03 

Heptane 

1.387 

0.0391 0.1331 

0.1621 

2.678 

1.036E02 

Tetradecane 

1.429 

0.0242 0.1184 

0.1766 

1.737 

2.573E01 


For the equilibrium case, Q - 0 and a - 1. Any values of 
a other than unity will signify deviation from the equilibrium 
situation that needs to be evaluated because of the sensitivity 
of the system to extremely small variations in the local tem- 
perature. The experimental section of the article and the 
subsequent discussion will present three experimental cases, 
which were very close to equilibrium. The results will demon- 
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strate the utility of the model and the evaluation of the val- 
ues of a. 

After multiplying both sides of Eq. 10 by Zdry'dZ and inte- 
grating (C, is the constant of integration). 



The boundary condition used for the completely wetting case 
is 


For r\ = a 


d ~T) di) 

~dZ^~°~dZ’ 


( 12 ) 


Note that this can be viewed as both an artificial boundary 
condition for a nonequilibrium system for which a > 1 and 
the actual boundary condition for an equilibrium system for 
which a =■ 1. 

Using the boundary condition, the slope of the meniscus 
can be expressed as 


dS yp I 2a* 8 


(13) 


Hence if the curvature at the thicker end of the meniscus, 
K„ along with fl, 5 0 , and a are known, the slope of the 
meniscus can be directly calculated as a function of the film 
thickness, using only the augmented Young- Laplace equa- 
tion. The minus sign in Eq. 13 is indicative of the fact that for 
the reference frame selected, the meniscus slope should al- 
ways be negative (film thickness decreases as distance in- 
creases). In a recent publication concerning a different exper- 
imental design, DasGupta et al. (1994) have demonstrated 
that the augmented Young-Lapiace equation can be success- 
fully used to evaluate the slope profiles for isothermal cases 
from known values of B, 5 n , and a. 

For our present study we obtain the curvature in the thicker 
portion of the liquid meniscus in the following way. As can 
be seen from Figures 6 and 7, the isothermal film profile in 
the thicker portion of the meniscus approximated a parabola 
(5 * ax 1 , with a being a constant and the curvature equal to 
2a). The values of 5 1/1 and the corresponding jc, for S > 0.5 
jim. were fitted to a straight line using linear regression. The 
result was an excellent fit in each case ( R 2 values higher than 
0.99) and from the regression routine the slope of the result- 
ing straight line (equal to a 1 '' 2 ) and K x (equal to 2a) were 
readily obtained. 

Herein we compare the slope of the meniscus obtained by 
numerical analysis of the experimental data for a system with 
unknown values of Hamaker constant (or dispersion con- 
stant, depending on the adsorbed film thickness) and the slope 
predicted by the augmented Young-Lapiace equation (Eq. 
13). The slope predicted by Eq. 13 is a function of a and the 
a corresponding to the closest match between these two 
slopes is selected to determine the value of the dispersion 
constant for the CVBT, in situ , for our isothermal experi- 
ments, as will be apparent in the next sections. Finally, we 
compare the experimentally obtained values of the dispersion 


constants and the disjoining pressures with the predictions 
from the DLP theory ( Dzvaloshinskii et al., 19M) for our ex- 
perimental system. 


Results and Discussion 

The objective of this work is to characterize the interfacial 
force field by evaluating the dispersion constant, in situ . for 
the solid-liquid-vapor system. Once the film thickness pro- 
file. including the adsorbed film thickness, <5 n , were obtained, 
the curvature at the thicker end of the meniscus. K. c . was 
determined. Equation 13 was then used to calculate the slope 
of the thickness profiles as a function of a and the a corre- 
sponding to the closest match between the experimental 
slope, obtained by numerical differentiation of a cubic spline 
fitted to the data, and Eq. 13 was used to determine the dis- 
persion constant for the specific liquid-solid system. 

Figures 9 and 10 show two pictures of the close match be- 
tween the experimental slope and the slope obtained by the 
solution of the augmented Young-Lapiace equation. The 
solid line represents the slope obtained by numerical differ- 
entiation of the data, whereas the triangles represent the so- 
lution of the augmented Young-Lapiace equation for the 
specific value of a at every tenth data point. In Table 2 the 
results for pentane, which are a strong function of the value 
of a, are summarized. Similar results were also obtained for 
heptane. 

The experimental values of the dispersion constants for the 
systems are presented in Table 3. The dispersion constants 
are calculated from the known values of a corresponding to 
the minimum error. The values of B for each liquid are rea- 
sonably close to each other. We also calculated the disper- 
sion force and the theoretical values of B from the DLP the- 
ory (Lifshitz, 1955; Dzyaloshinskii et al., 1961) for quartz- 
pentane-vapor and quartz-heptane-vapor systems. To calcu- 
late the dispersion force from the DLP theory, the same 
methodology presented by Truong and Wayner (1987) and 
Gee et al. (1989) is used. The dielectric functions of the liq- 
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Figure 9. Experimental slope ve. that obtained from the 
augmented Young-Lapiace equation (Eq. 13). 
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Figure 10. Experimental slope vs. that obtained from thn 
augmented Young-Laplace equation (Eq. 
13). 

uids and substrates are obtained from spectroscopic optical 
data given in Table 4 as described by Hough and White (1980). 
The results are presented in Table 3 and in Figure U. Figure 
1 1 is a plot of the disjoining pressure as a function of the film 
thickness for pentane and heptane wetting on quartz, and for 
tetradecane on SFL6. The solid line is from DLP theory, 
whereas the symbols are from the experiments. Since the val- 
ues of the dispersion forces for the two alkanes on quartz, as 
predicted by the DLP theory, are about the same, only the 
theoretical values for pentane are presented. Both Table 3 
and Figure 11 clearly demonstrate that our experimentally 
obtained values of B and the disjoining pressures for the two 
liquids in the quartz cell are close to those predicted from 
the exact DLP theory. Comparable results for tetradecane on 
SFL6 are presented in Figure 11 and Table 5. However, the 
DLP calculations for SFL6 are based on incomplete optical 



Figure 11. DLP theory for an ideal surface vs. experi- 
ments for the adsorption of pentane and 
heptane on quartz ( — ) and tetradecane on 
SFL6 based on UV data only ( ). 


data. We note that the DLP theory predicts the dispersion 
forces for an ideal system, whereas the surface energies of 
real surfaces can be quite different because of the presence 
of impurities or adsorbed water layers. The importance of an 
in situ force field characterization technique cannot therefore 
be overemphasized. 

Similar results were obtained by Gee et al. (1989) from 
adsorption studies of n-alkanes on quartz for high disjoining 
pressure regimes (5 <> 40 A) and by Blake (1975), who mea- 
sured the disjoining pressure as a function of film thickness 
for n-octane and n-decane on a-alumina. Blake’s techniques 
demonstrated that Lifshitz theory correctly predicts the re- 
sults in the low disjoining pressure regimes ( 5 was between 
250 A and 800 A), whereas our present study is concerned 
with adsorbed film thicknesses of about 200 A. For complete- 
ness, we also note that recent studies by Beaglehole et al. 


Table 2. Characteristics of Pentane Meniscus 




Angle of Inclination - 5.67* 


Angle of Inclination — 9.92° 



^0 

(nm) 

K x x 10* 1 
(ra* 1 ) 

B 

a (I'm) 

So 

(nm) 

K x x 10* J 
(m* ') 

a 

B 

(Jm) 

Point A 

22.0 

2.512 

1.12 1.43 X 10* 29 

19.0 

3.044 

1.04 

0.72x10* 29 

Point B 

20.0 

3.147 

1.04 0.91x10* 29 

20.0 

4.186 

1.12 

1.63x10 


Table 3. Experimental and Theoretically Calculated Values of B( 8) for Pentane and Heptane on Quartz 



n-pentane 




n-heptane 


Angle 

(°) 

«o 

(nm) 

^e»pt 

(Jra) 

&DLP 

(Jm) 

Angle 

(°) 

«o 

(nm) 

^expl 

(Jm) 

&DLr 

(Jm) 


22.0 

1.43 XlO* 29 

1.14X10* 29 


24.0 

2.60 xlO" 29 

1. 07 x 10* 29 

5.67 

20.0 

0.91 XlO* 29 

1.11 x 10* 29 

6.24 

23.0 

1.57 xlO* 29 

1.05 xlO* 29 


19.0 

0.72X10* 29 

1.09X10" 29 


23.0 

2.41 xlO* 29 

1.05 XlO* 29 

9.92 

20.0 

1.63 xlO* 29 

1.11 xlO* 29 

9.92 

22.0 

2.71 xlO* 29 

1.04 x 10* 29 
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Table 4. Optical Constants Used the Calculation of 
Dispersion Force from DLP Theory 


4 



C ut 

(o Ll 

( rad/s) 

C,„ 

10 IR 

( rad/s) 

Pentane 

0.819 

1.877 x lO 16 

0.025 

5.54 X 10 14 

Heptane 

0.898 

1.870 x iO 16 

0.025 

5.54 x 10 1J 

Tetradecane 

1.01 1 

1.846 x 10 16 

0.025 

5.54X10' 4 

Quartz 

1.359 

2.032 x lO 16 

1.93 

2.093 x I0 14 

SFL6 

rn 

0 

r i 

1.191 x 10"’ 




(From Hough and White. IdSO). 


(1991) and Beaglehole and Christenson (1992) demonstrate 
the limitations of the DLP theory in predicting the interac- 
tions in very thin liquid films (below 20 A), and their results 
suggest that structural effects are present in some adsorbed 
films at room temperature, provided the substrates are 
smooth and homogeneous. 

Effects of evaporation on a 

The augmented Young- Laplace equation can be written in 
the following way 


1 1 d 2 S 

A/Sq 17 3 at 3 K « dx 2 


(14) 


where 0 is the dimensionless pressure. Note that we are us- 
ing a different expression for the disjoining pressure (equal 
to A/Sq\ and hence the nondimensional augmented 
Young- Laplace equation is expressed in a slightly different 
way with B/S 0 substituted by A in the expression for a. Nev- 
ertheless, all the physical significances of a remain un- 
changed. This was done because we will be using some exper- 
imental data from another study and the magnitude of the 
adsorbed film thickness in that study warranted the use of 
the modified Hamaker constant instead of the dispersion 
constant. 

In the thicker portion of the meniscus (large tj), the previ- 
ous equation can be simplified to obtain 

a 3 *—, 7i-*®. (15) 


For an ideal equilibrium case, 0 — 1 and there is no evapora- 
tion, condensation and flow in the meniscus. Whereas, for a 
nonequilibrium case, 0 deviates from 1 and a larger than 
unity value of a is obtained when evaporation occurs. Since 
0 in Eq. 15 is equal to the ratio of the pressure jump across 
the liquid-vapor interface at large tj to the pressure jump at 


Table 5. Experimental and Theoretically Calculated Values 
of B( S ) for Tetradecane on SFL6 




Tetradecane 


Angle 

(°) 

So 

(nm) 

a 

B eX pt 

(J*m) 

Bdlp 

(J*m) 


27.0 

1.14 

5.87x10’ 29 

1.60x10- 29 

6.88 

27.0 

1.15 

7.12 x 10" 29 

1.60x10- 29 


EFFECTS OF EVAPORATION 

ON 

a 
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qmax (W/m 2 ) 

Figure 12. Effects of evaporation on the values of a for 
propanol and R113 menisci. 


rj - 1, we find that (1 - 0 tj) is a measure of the average di- 
mensionless shear stress at the liquid -solid interface due to 
liquid flow. The theoretical value of 0 for a given evapora- 
tive heat flux can be obtained using the equations presented 
by DasGupta et al. (1993b). Therefore, we can connect the 
very small value of a found herein with an extremely small 
amount of phase change. 

We reanalyzed the data for some of the nonisothermal (and 
a few isothermal) cases observed in a different experimental 
setup designed to study evaporation. In those experiments a 
capillary feeder was used to study and model evaporating 
menisci of various liquids on a silicon substrate (DasGupta et 
al., 1993b). The values of 6 0 were accurately measured using 
ellipsometry. A theoretical model based on the augmented 
Young-Laplace equation and modified Kelvin -Clapeyron 
equation was used. The numerical solution of the proposed 
model and comparison of the results with experimental data 
were used to characterize the system and to obtain the di- 
mensionless pressure difference (0) and maximum evapora- 
tive heat flux. The maximum heat flux data were presented in 
a related publication as a function of the interfacial tempera- 
ture difference (DasGupta et al., 1994). Herein we utilize 
these results to calculate the values of a to demonstrate the 
sensitivity of a to evaporation. 

Figure 12 is a plot of a as a function of maximum evapora- 
tive heat flux for nonequilibrium propanol and R113 menisci. 
The figure clearly demonstrates that a increases with an in- 
crease in evaporation (departure from equilibrium). The first 
points for both propanol and R113 data are interesting to 
note. For both of these cases the heater was off, but as the 
values of a show, the case for propanol was very close to 
isothermality (a - 1.037 with a very small value for the heat 
flux), whereas the case for R113 is farther from equilibrium 
(a - 2.05) due to uncontrolled environmental conditions. The 
respective temperature differences between the solid and va- 
por were found to be approximately 10 -5 and 10 -4 K. This is 
an example of the sensitivity of the whole process. We also 
note that a small change in the value of a can signify a large 
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change in the evaporative heat flux. A value of a greater 
than one denotes the presence of a curvature gradient in the 
flow field and shear stress due to the resulting flow. 

The most important point of this study is the demonstra- 
tion that the dispersion constants for a constrained vapor 
bubble thermosyphon can be evaluated in situ at the start of 
the experiments. This is a definite improvement over previ- 
ous studies of micro heat pipes where little attention was given 
to the values of dispersion constants that are a strong func- 
tion of the previous history of the substrate used in the ex- 
periments. These constants characterize the interfacial force 
field and are profoundly important in the basic understand- 
ing of the operation and performance of a micro heat pipe 
and other small thermosyphons (which includes the micro- 
layer in boiling). 

Conclusions 

• The use of IAI in conjunction with a CVBT in the earth’s 
gravitational field under isothermal conditions was demon- 
strated. 

• The accuracy of the results are a function of the difference 
between the refractive indices of the solid and the liquid. 

• Procedures to measure the dispersion constant in situ for 
the vapor-liquid-solid system were developed and good 
agreement between theoretical and experimental values of the 
dispersion constant was obtained. 

• The augmented Young-Laplace equation accurately mod- 
els the change in the effective pressure at the liquid -vapor 
interface. 

• The sensitivity of the process and the physical significance 
of the parameter a in relation to the closeness to equilibrium 
were demonstrated. 
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Notation 

x - distance, m 

Z - dimensionless distance, Eq. 7 
t; - dimensionless film thickness, 5/3 0 
s - solid 

Y- L - Young-Laplace equation 
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